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ABSTRACT 



Context. The contrast of granulation is an important quantity characterizing solar surface convection. 

Aims. We compare the intensity contrast at 630 nm, observed using the Spectro-Polarimeter (SP) aboard the Hinode 

satellite, with the 3D radiative MHD simulations of Vogler & Schiissler (2007). 

Methods. A synthetic image from the simulation is degraded using a theoretical point-spread function of the optical 
system, and by considering other important effects. 

Results. The telescope aperture and the obscuration by the secondary mirror and its attachment spider, reduce the 
simulated contrast from 14.4% to 8.5%. A slight effective defocus of the instrument brings the simulated contrast down 
to 7.5%, close to the observed value of 7.0%. 

Conclusions. A proper consideration of the effects of the optical system and a slight defocus, lead to sufficient degradation 
of the synthetic image from the MHD simulation, such that the contrast reaches almost the observed value. The 
remaining small discrepancy can be ascribed to stray light and slight imperfections of the instrument, which are difficult 
to model. Hence, Hinode SP data are consistent with a granulation contrast which is predicted by 3D radiation MHD 
simulations. 
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1. Introduction 

The root-mean-square of the normalized continuum inten- 
sity fluctuations, within an area on the solar disk, is deter- 
mined mainly by the intensity variation between the bright 
^ granules and the darker intergranular lanes, and thus is 
K> , usually referred to as the granulation contrast. It is a key 
^ ■ property of solar surface convection because it is connected 
d I with the temperature difference between rising (granules) 
and descending gas masses (intergranules), and thus related 
to the efficacy of the convective energy transport. 

Reliable measurements of the granulation contrast are 
notoriously difficult, since the observed contrast suffers sig- 
nificantly from image degradation by the optical system 
and, most importantly in the case of ground-based tele- 
scopes, by seeing and straylight effects due to the terrestrial 
atmosphere. To deduce the 'true granulation contrast' from 
measurements, we need to deconvolve the observed images 
by considering the modulation transfer through the optical 
setup (telescope, instruments, detector, etc.), in addition 
to through the terrestrial atmosphere. The quantitative ef- 
fects of the latter are poorly understood and it is therefore 
unsurprising that the reconstructed values of the granula- 
tion contrast, documented i n the literature, cover a broa d 
range (see, e.g., Table 2 in 'Sanchez Cuberes et al.]|2QQQl ). 
This remains true even if additional constraints, such as 
measurements of the intensity profile across the lunar limb 
during a partial eclipse, are used to derive the effective 



point-spread function (e.^., Lew Il971; [Peubner fc Mattid 
119751 : lNordlundl ll984': "Sanchez Cuberes et al."'2000). 

In the case of balloon-borne stratospheric observations, 
the influence of the atmosphere is negligible. Even under 
these conditions, however, t he scatter of the reconstructed 
contrast values is la r pe (iBahng fc Schwarzschildl 1961 



Pravdvuk et al. 1974; Altrock 1976; Edmonds 
1977; Witt mann fc Mehltretter 1977; Schmi dt et al 



linkl. 
[19? 



Wit t manni ll 98 if ) , presumably because the p roper consider- 
ation of the instrumental effects is nontrivial ([Durrant et al.l 
1983). 

Although the range of contrast values derived from 
observations is quite wide, the values predicted by 
3D radiativ e HP /MHD sini ulations are significantly 
hig her fe.g. IStein fc Nordlundl l2QQQl : I Vogler et all l2QQ5l : 
Wed emever-Bohm 2QQ7I ). This is unsurprising, since the 
horizontal resolution of such simulations (10-40 km) is 
far better than what can be achieved observationally at 
present. To compare the observed contrast values with 
predictions from simulations, the synthetic images derived 
from simulations must be degraded in the same way as 
the observations, taking into account the effects of the op- 
tical system and, in the case of ground-based telescopes, 
the terrestrial atmosphere. This procedure is prone to the 
same uncertainties as the deconvolution of observations; 
the result depends strongly on assumptio ns about the na- 
ture of atmospheric seeing and straylig ht ([Nordlundl 1 19841 : 
ISchiissler et~al]l2QQ3l : [Rvbak et al.ll2QQ6h and we are unable 
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to infer if the simulation predictions are consistent with ob- 
servations, or otherwise. However, the good agreement be- 
tween observed and simulated spatially averaged spectral 
line bisectors indicates that the simulated intensity con- 
trast s are probably close to actual values (|Asplund et alJ 

A new era has begun with the launch of the 50-cm 
Solar Optical Teles cope (SOT, Tsuneta et al. I j2008) on 
the Hinode satellite (jKosugi et al.ll2QQ7l ). The good perfor- 
man ce and low stray lig:ht level of the spectro-polarimeter 
(SP. lLites et al.|[2QQll )and complete absence of atmospheric 
effects, enable a far more reliable determination of the gran- 
ulation contrast at 630 nm, the wavelength at which the SP 
operates. In this Letter^ we compare the intensity contrast 
of a Hinode SP continuum map of the quiet Sun with pre- 
dictions of MHD simulation. 



2. Observations and simulation data 

We determine the observed contrast from a map of the con- 
tinuum intensity at 630 nm wavelength obtained on Jan. 16, 
2007 (12:10:10 - 13:36:49 UT), using the scan mode of the 
Hinode SP with an exposure time of 4.8 s per slit position. 
The map covers a field of 163'' x 164'' of 'quiet' Sun, close 
to disk center. The slit width, the sampling step size, and 
the CCD pixel pitch all correspond to 0.16". The data were 
reduced using stan dard routines that correct for various in- 
strumental effects (|Lites et al.l 120071 ). After reduction, we 
find an rms contrast of the continuum intensity of 7.0%. 

The 3D MHD simulation snapshot that we consider rep- 
resents the satura ted (statistically sta t ionar y) state of the 
dynamo run C of lVogler &: Schiissled (|2007r ). The compu- 
tational box contains 648 x 648 x 140 cells; it corresponds 
to a physical size on the Sun of 4.86 x 4.86 Mm^ in the 
horizontal and 1.4 Mm in the vertical direction, the latter 
ranging from about 900 km below to 500 km above contin- 
uum optical depth unity at 630 nm (reso = !)• The simu- 
lation has been run with non-grey radiative transfer using 
the MURaM code (Vogler 2003; Vogler et al. 2005). With 
a horizontal grid-cell size of 7.5 km, this simulation had one 
of the highest resolutions achieved so far. The snapshot that 
we use, has an average unsigned vertical magnetic field of 
about 7 G at reso =0.1, which we consider as a reason- 
able representation of the 'quiet' Sun. The corresponding 
original (unsmeared) continuum image at 630 nm shows an 
rms contrast (standard deviation divided by mean value) 
of 14.4%. The particular choice of snapshot was not criti- 
cal for the results presented here. Considering a number of 
other snapshots and simulations of different spatial resolu- 
tion or amount of magnetic flux in the simulation box, we 
found contrast values that were typically between 14% and 
15%. 



3. Modeling of the system PSF 

The difference between the observed and original contrast 
of the simulation amounts to more than a factor of two 
(7.0% vs. 14.4%). Since there is no atmospheric distortion 
of the image, we investigate whether instrumental effects 
alone can account for this significant degradation. 

The intrinsic resolution of the simulated images consid- 
ered here (equivalent to 0.01" pixel size) is considerably bet- 
ter than the resolution of the 50-cm Hinode SOT; the degra- 
dation of the synthetic image based on the simulation result 
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Fig. 1. Change of the Modulation Transfer Function when 
the effects of the different parts of the optical system are 
sequentially taken into account: telescope aperture ('Airy'), 
central obscuration, secondary mirror and spider, sampling 
(CCD), and a defocus of 1.5 mm (about 9 steps of the focus 
mechanism) . 



can therefore be modeled by applying an effective point- 
spread function (PSF). We convolved our simulated inten- 
sity maps in the 2D spatial domain with different PSFs, 
which represented a more and more realistic optical sys- 
tems. To this end, we modeled the telescope and the trans- 
fer optical path down to the spectrograph entrance slit us- 
ing the commercial optical design software ZEMA^xEl- While 
the telescope was modeled with the nominal SOT surface 
parameters, the transfer (refocussing) path was modeled in 
paraxial approximation, which, however, has no influence 
on the final results of our calculations. ZEMAX calculates 
the PSF as the Fourier transform of the wavefront in the 
exit pupil of the system. 

We describe the contrast degradation by the optical sys- 
tem using the Modulation Transfer Function (MTF), which 
is the Fourier transform of the point-spread function. In this 
way, the various effects can be considered step by step, as 
indicated below. Figure [T] shows the corresponding change 
of the system MTF. 

Telecope aperture: The Hinode SOT is an //9.1 Gregory 
system with an aperture of 50 cm, which corresponds to a 
cut-off at 1/0.26". Higher spatial wave numbers are not 
transmitted by the system, while low to intermediate wave 
numbers are transmitted with relatively high contrast. 

Central obscuration and spider: Important contributors 
to the contrast reduction are the central obstruction of the 
telescope by the secondary mirror, which produces a linear 
obscuration of 0.344 = 17.2 cm/50 cm, and the diffraction 
by the three spider elements, which are of 4 cm width each, 
holding the secondary. The effect of the central obstruction 
is that the highest spatial wave numbers are largely un- 
affected, while the intermediate wave numbers are signifi- 
cantly damped. The spider structure adds to this damping. 
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Korrralized intensity 

Fig. 2. Normalized histograms of the continuum-intensity 
values at 630 nm, obtained for 35 equal bins between 0.6 
and 1.4 of the mean intensity. Shown are results for the orig- 
inal synthetic image from the simulation snapshot (outer- 
most, solid curve) and for the degraded images correspond- 
ing to the MTFs shown in Fig. [1] (same line styles as used 
there). The two thick inner curves represent the degraded 
synthetic image assuming a defocus of 1.5 mm (thick dashed 
curve) and the continuum map observed with Hinode SP 
(thick solid curve). 



Spatial sampling: Since the data considered here are 
obtained using the spectro-polarimeter, we have to con- 
sider the sampling and integration effects of the slit and 
the CCD detector. The sampling and the integration ef- 
fects of the detector are approximated by powers of sine 
functions (Boreman 2001). 

Defocus: Low-order optical aberrations damp the MTF 
at low spatial wave numbers. Since the Hinode Focal Plane 
Package (FPP) is a complex system, it is out of the scope of 
this Letter to model in detail all relevant contributions from 
optical, optomechanical, and electronic parts of the system. 
We therefore consider only the effect of a defocus in the 
plane of the SP slit. It turns out that a small amount of de- 
focus can reduce significantly the intensity contrast, while 
leaving the fine-scale resolution of the system unaffected. 
The focus of the Hinode science instruments is controlled 
by means of a common reimaging lens that can be shifted 
in a range of ±25 mm, with a step of 0.17 mm. 

In addition, we investigated the effect of the (lossy) 
JPEG compression which was applied to the images. The 
change of the rms contrast turned out to be negligible, even 
for high compression factors. 

4. Degradation of the simulation data 

To compare with the Hinode result, we convolved the inten- 
sity maps calculated from the simulation with the different 
PSFs resulting from modeling the Hinode SOT as described 
in the previous section. After convolution, the degraded im- 
ages were binned to represent the sampling and integration 
effects of the finite slit width, step size of the spatial sam- 
pling, and the pixel size. 



Table 1 . Values of the rms intensity contrast after applying 
the PSFs corresponding to the MTFs shown in Fig. [TJ 



effects taken into account rms [%\ 



none 


14.4 


Airy 


10.9 


-h central obscuration 


9.6 


-h spider 


8.7 


-h CCD 


8.5 


-h defocus (1.5 mm) 


7.5 



Figure [2] shows the evolution of the histogram of inten- 
sity values in the simulation image as the various effects of 
the optical system are introduced. The rms contrast val- 
ues for the various steps are listed in Table [TJ The basic 
telescope effects (primary aperture, obscuration, and spi- 
der) reduce the rms contrast to 8.7%, presumably since the 
small-scale intensity structure in the intergranular lanes is 
lost. There is almost no further contrast degradation due 
to sampling and rebinning (CCD), which indicates that the 
remaining contrast is determined mainly by spatial scales 
that are significantly larger than the critical sampling limit 
of the diffraction pattern. Once the fine detail in the gran- 
ular lanes is no longer resolved, the contrast is dominated 
by the intensity differences between the granules and the 
average intensity of the intergranular lanes, which is almost 
unaffected by the integrating effect of the detector pixels. 
Similar valu es of the contrast reduction up to this step were 
obtained by lOrozco Suarez etld] (|2007[ ) . 

The difference between the rms contrast degraded so far 
(8.5%) and the observed contrast (7.0%) is still significant. 
Another factor that contributes to contrast degradation is 
a slight defocus of the SP, which affects mainly the low to 
intermediate spatial wave numbers. We have no direct infor- 
mation about the amount of defocus for the dataset consid- 
ered here; we, however, obtained empirical evidence about 
the effect of defocus on the intensity contrast using a num- 
ber of SP datasets for different focus positions, taken only 
an hour earlier than the map analyzed here. Apart from 
the exposure time of 1.6 s per slit position and the smaller 
size of the maps, all other observational parameters are the 
same. For every focus position, a map that covers a field of 
3.2^' X 82^' was generated. The corresponding contrast val- 
ues given in Fig. [3] confirm the rather sensitive dependence 
on focus position. The figure also shows the dependence of 
the contrast of the simulation image on the defocus value. 
The large observational map with a contrast of 7.0% (indi- 
cated by the asterisk) was taken with a defocus of 8-9 steps 
(1.36 mm - 1.53 mm) from the optimal position in Fig. [3l 
The degraded simulation image corresponding to about this 
defocus value (1.5 mm) has a contrast of 7.5%. 

As an illustration of the image degradation, Fig.[4]shows 
the original continuum image from the simulation (left 
panel, the periodic simulation box is drawn fourfold for 
the sake of better visibility) and the final degraded image 
(middle). A subset of the observational map of the same 
size is shown in the right panel, using the same gray scale. 

The remaining deviation of 0.5% from the observed 
value of about 7%, and the general shift between the curves 
for simulation and observation in Fig. [3l can be ascribed 
to various degrading factors that are not included in our 
analysis. These include the effects of stray light, of all other 
low-order optical aberrations apart from defocus, and of the 
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Fig. 4. Continuum images at 630 nm from the simulation snapshot with original resolution (left) and after degrading 
(middle)^ in comparison with a detail of the observed Hinode/SP map of the same size (right). The periodic simulation 
box is shown fourfold for better visibility. 
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Fig. 3. Intensity contrast as a function of telescope defo- 
cus. Shown are results for a SP data set used taken on Jan. 
16, 2007, 10:54:14-11:06:52 UT (diamonds) and the con- 
trast values for the degraded simulation image (crosses). 
The value for the large map used to compare with the sim- 
ulation is indicated by the asterisk, assuming that the lo- 
cation of the optimal focus was the same as for the other 
dataset. This also gives an indication of the observational 
uncertainty. 

read-out electronics (such as transfer efficiency and pixel- 
to-pixel crosstalk), which reduce the contrast, but are diffi- 
cult to model reliably. To illustrate the extreme sensitivity 
of the granulation contrast to low-order aberrations and 
straylight we calculated, as an example, in ZEMAX the 
effect of astigmatism and coma, which were equally dis- 
tributed in such a way that the total system shows an rms 
wave- front error of 0.044 wavelengths. This corresponds to 
the number given by Suematsu et al. (2008) and represents 
a system with a remarkable Strehl ratio of 0.93. This small 
effect is nevertheless sufficient to reduce further the con- 
trast from 8.5% to 8.1% in optimal focus position. If no 
further degrading effects are present, a straylight level of 



only 4.7% would be sufficient to bring the contrast down to 
the observed value of 7.8% (for optimal focus), and to 7.0% 
for the presumptive defocus value of the dataset studied 
here, respectively. 

5. Conclusion 

We find that consideration of the basic optical properties of 
the Hinode SOT/SP system and a slight defocus are suffi- 
cient to bring the degraded contrast of a 3D radiative MHD 
simulation and the observed rms continuum contrast at 
630 nm, almost into agreement. The remaining discrepancy 
can be ascribed to the combined, minor effects of stray- 
light, other low-order optical aberrations, and instrument 
electronics. Hence, Hinode SP data are consistent with a 
granulation contrast at 630 nm of 14-15%, as predicted by 
the simulations. 
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